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ABSTRACT 



Aims. We present a study of the orbit of the pre-main-sequence binary system GG Tau A and its relation to its circumbinary disk, in 
order to find an explanation for the sharp inner edge of the disk. 

Methods. Three new relative astrometric positions of the binary were obtained with NACO at the VLT. We combine them with data 
from the literature and fit orbit models to the dataset. 

Results. We find that an orbit coplanar with the disk and compatible with the astrometric data is too small to explain the inner gap of 
the disk. On the other hand, orbits large enough to cause the gap are tilted with respect to the disk. If the disk gap is indeed caused 
by the stellar companion, then the most likely explanation is a combination of underestimated astrometric errors and a misalignment 
between the planes of the disk and the orbit. 
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1. Introduction 

GG Tau is a young q uadruple system consisting of two binaries 
dLeinert et al.L 1 19931) . GG Tau A is a pair of low-mass stars sep- 
arated by about 0.25". GG Tau B, located 10.1" to the south, 
is wider (1.48") and less massive. A circumbinary disk around 
GG Tau A has been extensively studied. It was spatially re- 
solved in both near infrared and millimeter wavelength domains. 
A detailed analysis of the velocity maps of the disk found that 
it is in Kepleria n rotation and constrai ned the central mass to 
1.28 + 0.07 Mm (iGuilloteau et all [19991 

So far, orbital motion has not been detected in the GG Tau B 
binary because of its long period. However, the relative motion 
of the components of GG Tau A has been observed for several 
years and has a lready resulted in several orbit determinations 
(iMcCabe et all 1200a iTamazian et ail l2002t iBeust & Dutrevl 
120051) . Because only a limited section of the orbit has been 
observed, the authors generally have assumed that the orbit is 
coplanar with the circumstellar disk (i = 37° + 1°). The result- 
ing orbital parameters were all quite similar to each other, with 
a semi-major axis of about 35 AU. 

The presence of the binary would be an obvious explanation 
for the rather sharp inner edge of the disk located at 180AU. 
The ratio of the inner radiu s of the disk and the semi- major 
axis is about five. However, Artvmowicz & Lubowl ([1994) stud- 
ied the effect of binary systems on their circumbinary disks and 
found that this ratio should range from about 1.7 (for circular 
orbits) to about 3.3 (for hig hly eccentric binaries, e = 0.75). 
Beus t & Dutrevl (120051 [2006) carried out a similar study specifi- 
cally for GG Tau A and came to the same conclusion. The binary 
orbit cannot explain the gap in the circumbinary disk, unless its 



* Based on observations collected at the European Southern 
Observatory, Chile, proposals number 072.C-0022, 078.C-0386, and 
384.C-0870 



semi-major axis is about twice as large as indicated by the astro- 
metric data available. 

In this paper, we present new relative astrometric measure- 
ments of GG Tau A and derive estimates for its orbital parame- 
ters, with and without the assumption that binary orbit and cir- 
cumbinary disk are coplanar. 



2. Observations and data reduction 

Astrometric measurements of GG Tau A have been pub 
lished by se v eral authors dDuchene et all 2004t Ghez et al 



19951 [1997 



Leinert et al 



lHart igan &Kenvon, 2003; Kristetal 
19931: McCabe et all l2002t iRoddier etal 
Tamazian et all 120021; IWhite & GhezT 



2002 



1996 



2001; Woitas et al 



20011) . see IBeust & Dutrevl (120051) for an overview. Here we 
report on new observations obtained with NAOS/CONICA 
(NACO for short), the adaptive optics, near-infrared camera 
at the ESQ Very Large Telescope onCerro Paranal, Chile 
dRousset et all 12003b iLenzen et all l2003h . GG Tau was ob- 
served on December 13, 2003 (PI: Leinert), November 20, 2006 
(PI: Ratzka), and October 5, 2009 (PI: Kohler). We use only 
imaging observations in the K„ photometric band for the orbit 
determination. Integration times were 85 sec per image in 2003, 
24 sec in 2006, and 60 sec in 2009. In 2003 and 2006, we took 
four images with the star at different positions on the detector 
to facilitate creation of a median sky image. In 2009, 12 images 
were recorded at four positions. 

The NACO images were sky subtracted with a median sky 
image, and bad pixels were replaced by the median of the clos- 
est good neighbors. Finally, the images were visually inspected 
for any artifacts or residuals. Figure Q~]shows an example of the 
results. 

The Starfinder program dDiolaiti et all 120001) was used to 
measure the positions of the stars. The positions in several im- 
ages taken during one observation were averaged, and their stan- 
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Table 1. Astrometric measurements of GG Tau Aa - Ab 



Fig. 1. Image of GG Tau A obtained with NACO in October 
2009. The separation between the two components is about 
250 mas. 



dard deviation used to estimate the errors. To derive the exact 
pixel scale and orientation of the detector, we took images of 
fields in the Orion Trapezium during each observing campaign, 
and reduced them in the same way as the images of GG Tau. The 
measured positions of the cluster stars were compar ed with the 
coordinates given in McCaughrean & Staufferl(fl99l . The mean 
pixel scale and orientation were computed from a global fit of all 
star positions. The scatter of values derived from star pairs was 
used to estimate the errors. The errors of the calibration are usu- 
ally comparable to or larger than the errors of the measured posi- 
tions of the science target, indicating the importance of a proper 
astrometric calibration. 

The calibrated separations and position angles of GG Tau A 
appear in Table Q] together with the data taken from the litera- 
ture. If one or both components of GG Tau B were within the 
field-of-view, then we also measured their positions. The results 
appear in Table [2] The main conclusion is that there has been no 
significant change i n the relative positi on since the first measure- 
ment published in lLeinert et al.l dl993b . 



3. Determination of orbital elements 



iMcCabe et ail d2002l) and iBeust & Dutrevl d2005l) have deter- 
mined the orbital elements of GG Tau Aa-Ab from the aver- 
age position a nd velocity of the com panion. Together with the 
system mass dGuilloteau et al.L Il999h . position and velocity on 
the sky comprise five measurements. Since orbital elements are 
seven unknowns, their computation requires the additional as- 
sumption that the orbit and the circumbinary disk are coplanar. 

In this work, we employed a different approach. We fit or- 
bit models to the observations and searched for the model with 
the minimum x 2 - in t he end, we wante d to use a Levenberg- 
Marquardt algorithm dPress et al.l 1 1992b . However, the results 
of this algorithm depend strongly on the chosen start values. To 
avoid any bias for a particular orbit, we carried out a preliminary 
fit that consists of a grid search in eccentricity e, period P, and 
time of periastron To- Singular value decomposition was used to 
solve for the remaining four elements. The result is a grid of x 2 
as function of e and P. Since we were interested in the semi- 
major axis a of the orbit, this was converted onto a a-e-grid by 



Date (UT) 



d [mas] 



PA[°] 



Reference 



1990 Nov 2 255. + 10. 9. ±2. Leinert et al. (1993) 

19910ct21 260. + 10. 2. ±1. Ghez et al. (1995) 

1993 Dec 26 260. ±10. 3. ±2. Roddie r et al. (1996) 

1994 Jan 27 246. ±4. 357.8 ± 0.4 Woitas et al. (2001) 
1994 Jul 25 250.2 + 2.6 358.8 ± 0.45 Ghez et al. (1997) 
1994 Sep 24 258. ±4. 357. ±2. Ghez et al. (1995) 
1994 Oct 18 242. ±3. 0.9 ±0.5 Ghez et al. (1995) 

1994 Dec 22 239. ±5. 357.2 ± 2. Roddie r et al. (1996) 

1995 Oct 8 247. ±4. 356.9 ± 0.7 Woitas et al. (2001) 

1996 Sep 29 245. ±4. 355.5 ± 0.4 Woitas et al. (2001) 

1996 Dec 6 243.6 ± 4.6 354.9 ± 1.3 White & Ghez (2001) 

1997 Sep 27 250. ±3. 354.3 ± 1. Krist et al. (2002) 
1997 Oct 10 248. ±2. 353.9 ± 0.4 McCabe et al. (2002) 

1997 Nov 16 247. ±5. 353.6 ± 0.4 Woitas et al. (2001) 

1998 Oct 10 260. ±4. 350.7 ± 0.4 Woitas et al. (2001) 

2001 Jan 21 248. ± 14. 348.6 ± 2.4 Hartigan & Kenvon (2003) 

2001 Feb 9 245. ±4. 348.7 ± 0.3 Tamazian et al. (2002) 

2002 Dec 12 250.7 ± 1.5 346.0 ± 1.5 Duchene et al. (2004) 

2003 Dec 13 250.7 ± 0.8 344.2 ±0.1 this work 
2006 Nov 20 252.3 ±0.7 339.0 ±0.1 this work 
2009 Oct 5 252.5 ±0.3 334.5 ±0.1 this work 



Table 2. Astrometric measurements of GG Tau B 



Date (UT) 


Pair 


d [arcsec] 


PA 


n 


2006 Nov 20 


Bb- 


-Ba 


1.460 


± 0.002 


134.9 


±0.1 




Aa- 


-Ba 


10.07 


±0.01 


185.4 


±0.1 


2009 Oct 5 


Aa- 


-Ba 


10.09 


±0.01 


185.5 


±0.1 



finding the orbit model with the closest a for each grid point. 
The grid spans a range from 20 to 200 AU in a, and from to 
0.99 in e. 

To convert the measured sep arations into AU, a distance of 
140 pc was adopted dEiiasLll978l) . 



3. 1 . Orbits coplanar with the disk 

First, we searched for an orbit matching all the information avail- 
able, i.e. the astrometric position, the total mass, and the ori- 
entation of the disk plane. We assumed that disk and orbit are 
coplanar, orbits without this constraint are discussed in the next 
section. 

The x 2 that we try to minimize is 

2 _ ST 1 ( r '.Obs _ r i.model \ 2 / M est - M mo d e l \ 2 



'disk 'model 



A/di 



isk 



disk 



model 



AQ, 



disk 



(i) 



where r l o b s and r,; mo d e i are the measured and predicted posi- 
tion at the time of observation i, and Ar, b s is the error of the 
measure ment. Here, M RSf i s the m easured system mass (1.28 + 
0.07 M Q , Guillotea u et all Il999t) . and M mode i the system mass 
predicted by the orbit model. Then, and Q d i s k are the incli- 
nation and position angle (PA) of the ascending node of the orbit 
of a disk particle, A/ d ; s k and AQdisk are their errors. The inclina- 
tion of the disk is 37 + 1° (Guilloteau et all Il999h . but it is in 
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Fig. 2-x 1 as function of a and e for orbit models that are coplanar 
with the circumbinary disk (Sect. IXTT l. The cross at a — 34 AU, 
e = 0.28 marks the minimum, the contour line around it encircles 
the 99.7 % confidence region (corresponding to 3<x in the case 
of normally distributed errors). The areas in various shades of 
gray are within the 5cr confidence region, i.e. orbit models in the 
white area can be excluded with 5cr confidence. 




semi-major axis [AU] 

Fig. 4. x 2 as function of a and e for orbit models that are not 
necessarily coplanar with the circumbinary disk (Sect. 13.2b . As 
in Fig. [2] the minimum^- 2 is marked by the cross, but the con- 
tour line around it encircles the 68.3 % confidence region (lcr). 
The areas in shades of gray are within the 5cr confidence region. 
The jagged shape of the contour line is most likely caused by 
numerical effects. 
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Fig. 3. Best-fitting orbit model if the orbit is constrained to be 
coplanar with the circumbinary disk (Sect. I3.lt . The observed 
positions are marked by their error ellipses and lines connecting 
the observed and calculated position at the time of the observa- 
tions. The new observations with NACO are marked by crosses. 
Their errors are too small to be discernible. The dash-dotted line 
indicates the line of nodes, the dashed line the periastron, and 
the arrow shows the direction of the orbital motion. 



Fig. 5. Three exemplary orbit models that fit the astrometric data 
and the system mass, but are not coplanar with the disk. The 
semi-major axes of the orbits are 35 AU, 85 AU, and 137AU. 
The last orbit has the minimal^ 2 , and its line of nodes is marked 
by the dash-dotted line and its periastron by the dashed line. The 
observed positions are marked by their error ellipses and lines 
connecting the observed and calculated position at the time of 
the observations. The new observations with NACO are marked 
by crosses. Their errors are too small to be discernible. 



retrograde rotation, so = 18 - 37 = 143°. The PA of the 
minor axis of the disk is 7 + 2° (Guilloteau et all fl999t) . there- 
fore Qdisk = 277 + 2° (the ascending node is defined as the point 
in the orbit where the objec t is receding from the observer most 
ranidlv. e.g. lHilditchll2001h . 



Equation[T]was minimized by a Levenberg-Marquardt algo- 
rithm dPress et all |T992). The starting points for the algorithm 
were taken from the preliminary fit described in the previous sec- 
tion. We kept a and e fixed to preserve the grid in these two vari- 
ables. The resulting x 2 distribution is depicted in Fig. [2] There 
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is a clear minimum at a — 34 AU and e = 0.28, while orbits 
with a > 36 AU can be excluded on the 3cr level. This is in per- 
fect agreement with previous orbit determinations. The reduced 
X 1 at the minimum is 3.05, which indicates a less-than-perfect 
fit. Figure [3] shows the orbit with the minimum^ 2 , together with 
the measurements of the relative positions, and Table |3]lists the 
orbital elements. 

To test whether the astrometric errors were underestimated, 
we repeated the procedure, but enlarged the errors of the ob- 
servations by a factor of 3. This lowers x 2 m general, but does 
not result in significant changes of the shape of the x 2 -plane as 
function of a and e. The best-fitting orbit has now a = 40 AU 
and e = 0.13. Also, because of the lower x 2 , many orbits with 
a > 36 AU (up to the end of the grid at a = 200 AU) are within 
the 99.7 % confidence region (which corresponds to 3<x in the 
case of a normal distribution). However, it appears unlikely that 
the authors of all astrometric data underestimated their errors by 
such a large factor, and orbits large enough to cause the disk gap 
are still only marginally consistent with the data. 

3.2. Orbits with no constraint on their orientation 

In this section, we remove the constraint that the orbit has to be 
in the same plane as the circumbinary disk. The only constraints 
are therefore the astrometric measurements, and the total mass 
of the binary. Then, x 2 is given by (using the same symbols as in 
Eq.EQ 



X 



-z 



'"/.obs '"/,model 



M est - M mode l \ 

AM est I 



(2) 



Figure [4] shows the result of minimizing the x 2 given by 
Eq.|2] The formal minimum is at a = 137 AU, e = 0.75 (TableO, 
with a reduced^ 2 of 3.2. It is highly unlikely that the true orbit 
has such a large semi-major axis and high eccentricity. However, 
the minimum is very shallow, and no semi-major axis larger than 
about 30 AU can be excluded, not even at the Icr level. Figure [5] 
shows the orbit with the formally minimal^ 2 , and two orbits that 
result in the best fit if the semi-major axis is held fixed at 35 AU 
and 85 AU, respectively. All three orbits fit the measured data 
reasonably well, demonstrating that the semi-major axis is not 
well constrained by the astrometric data. 



4. Discussion and conclusions 

If we require the orbit model to lie in the same plane as the cir- 
cumbinary disk, then orbits consistent with the astrometric data 
are not large enough to explain the gap in the disk. On the other 
hand, if we consider orbits that are not coplanar with the disk, 
then the astrometric data only provides a very weak constraint 
for the semi-major axis. This means that we can easily find or- 
bits that are consistent with the measured positions and with the 
size of the gap in the circumbinary disk. 

According to Artvmowicz & Lubowl dj994), an orbit with 
eccentricity e * 0.4 . . . 0.5 can open a disk gap that is a factor 
of about 3 larger than its semi-major axis. For our disk with an 
inner edge at about 180 AU, a semi-major axis of 60 AU would 
suffice. Figure [5] shows that orbits with a = 60 AU should have 
an eccentricity of 0.4 . . . 0.45 to match the astrometric data. We 
consider this to be the most plausible orbit, given the constraints 
from the astrometric data and the size of the disk gap. Its orbital 
elements appear in the rightmost column of Table [3] 



semi-major axis [AU] 

Fig. 6. Angle between the plane of the binary orbit model and 
the plane of the circumbinary disk. Plotted are the angles for 
all orbits within Icr of the best fit, i.e. all the orbits within the 
white contour line in Fig. |4] These are orbit models that were 
not constrained by the orientation of the disk. 



iBeust & Dutrevl (120051) have already discussed noncopla- 
nar solutions for the binary orbitQ Assuming a — 62 AU and 
e = 0.35, they find i = 125.4° and four solutions for Q: 114.3°, 
-65.7°, 74.5°, and -105.6°. The first solution differs from our 
most plausible orbit by 7° in inclination and 16.7° in Q, which 
is a reasonable agreement given the large uncertainties. 

How large is the misalignment between the plane of the orbit 
and the plane of the disk? Figure|6]shows the angle between orbit 
and disk plane as a function of the semi-major axis of the orbit. 
The relative astrometry of the two stars contains no information 
about the sign of the inclination (whether the orbit is tilted to- 
wards or away from the observer). In the creation of Fig. |6] we 
adopted in each case the inclination that resulted in the smaller 
angle between disk and orbit. 

Most orbit models are tilted by less than about 35° with 
respect to the disk. Our most plausi ble orbit is incli n ed by 
about 25°, a significant misalignment. IBeust & Dutrevl (2005) 
point out that the disk should show a warped structure if it 
is not coplanar with the binary orbit. This has not been de- 
tected, making such an orbit un likely, although it cannot be 
ruled out. Beust & Dutrev (2006) carried out some simulations 
of the dynamical b ehavior of the dis k for the noncoplanar or- 
bits found by Beus t & Dutrevl (120051) and a few possible orbits 
for the outer companion GG Tau B. In all the simulations with 
the orbit similar to ou r most plausible orbit (designated AA5 in 
IBeust & Dutrevl 120061) . the disk tends to assume an open-cone 
shape with an opening angle of ~ 30°. This state is reached after 
15 million years at the end of simulations. These results sug- 
gest that the GG Tau system we observe today is only a transient 

feature. 

GG Tau is only about 1 million years old dWhite & GhezL 

2001), so it is possible that we happen to observe the disk just 
before it dissolves into an open cone. We can only speculate how 
the system got into this unstable state. It is well known that stars 
experience strong gravitational interactions early in their life- 



1 They use a slightly different notation, where the position angle of 
the ascending node is replaced by the position angle of the projection 
of the rotation axis of the orbit onto the plane of the sky. The difference 
between the two position angles is exactly 90°. 
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Table 3. Parameters of the best orbital solutions. 



Orbital Element 


Orbit coplanar 


Orbit not 


most plausible orbit 




with disk 


coplanar w. disk 


tsee oect.[ftj 


Date of periastron To 


2477680 +2?o 


2460050 + _™ 


2463400 + _\™ 




(July 2071) 


(April 2023) 


(June 2032) 


Period P (years) 


162 -15 


i At\r\ + 1 7700 

1 400 

"""-1300 


403 _g 


Semi-major axis a (mas) 


243 -io 


977 


429 


Semi-major axis a (AU) 




137+17 

LJI -16 


60 


Eccentricity e 


0-28^ 


0.75 «j| 


44+0.02 

-0.03 


Argument of periastron a> (°) 


91 +4 

yl -13 


8 3 




PA. of ascending node CI (°) 


977 +2.0 
z ' ' -2.0 


318^° 


131!" 


Inclination i (°) 


143 +13 


128 « 


132.5 + _ l 2 ° 


Angle between orbit and disk 


0.02+1.9 


31.8 ±1.5 


24.9 ±1.7 



times, ev en if they form in small en sembles of only three to five 
stars (e.g. Ster zik & Durisenl Il998l) . These interactions can lead 
to catastrophic changes in binary orbits and even to the ejection 
of stars. The four stars in the GG Tau system would be enough 
to cause such events, unless they are in a stable configuration. 
Unfortunately, we have no kinematic information about the orbit 
of GG Tau B, which is not surprising, since we expect an orbital 
period on the order of 40000 years (based on the projected sep- 
aration of 1400 AU). It is conceivable that GG Tau has recently 
suffered a gravitational interaction and is currently in a transient, 
unstable state. However, a gravitational interaction that changes 
the orbit of GG Tau A should also have an effect on the circumbi- 
nary disk, making it highly unlikely that the disk could maintain 
the planar structure we see. 

On the other hand, the orbital elements derived from the 
astrometric data have rather large uncertainties. For example, 
the lcr confidence interval for the inclination of the orbit with 
a « 85 AU ranges from 115° to 158° (based on^ 2 as function 
of inclination). The errors of the angle between orbit and disk 
should be comparable, although not identical, since the angle 
between orbit and disk also depends on the orientation of the 
line of nodes. 

In summary, we do not have the final answer about the rela- 
tion between the orbit of GG Tau A and its circumbinary disk. 
An orbit coplanar with the disk could only cause the inner gap of 
the disk if the errors of the astrometric measurements are much 
larger than estimated. An orbit inclined to the plane of the disk 
would be compatible with both the astrometric data and the disk 
gap, but it should cause visible distortions in the disk structure. 
An explanation for the fact that no distortions in the disk have 
been detected could be that the orbit GG Tau A has only been 
changed recently, although any effect that can change the orbit 
of the stars should also disturb the structure of the disk. On the 
other hand, we should not forget the possibility that the gap in 
the disk is not related to GG Tau Ab, but some hitherto unknown 
companion. However, another companion would be pure specu- 
lation. 

The most likely explanation seems to be a combination of 
slightly underestimated astrometric errors and a (small) mis- 
alignment between the planes of the orbit and the circumbinary 
disk. More observations over a larger section of the binary orbit 
are needed. 

Acknowledgements. I thank the referee Herve Beust for his comments and sug- 
gestions that helped to improve the paper. 
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